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Description 

This invention relates to apparatus for determining the position and orientation of a remote object 
relative to a coordinate frame; and, more particularly, to radiating an electromagnetic field from the 

" 5 reference point, detecting the field at the remote object and analyzing the detected field to determine the 
position and, preferably, the orientation of the remote object. 

The use of orthogonal coils for generating and sensing magnetic fields is well known. For example, 
. such apparatus has received wide attention in the area of mapping magnetic fields to provide a better 
understanding of their characteristics. It has also been perceived that, if a magnetic field around generating 
10 coils can be very accurately mapped through the use of sensing coils, it might be possible to determine the 
location of the sensing coils relative to the generating coils based on what is sensed. However, a problem 
associated with doing this is that there is more than one location and/or orientation within a usual magnetic 
dipole field that will provide the same characteristic sensing signals in a sensing coil. In order to use a 
magnetic field for this purpose, additional information must therefore be provided. 

is One approach to provide the additional information required for this purpose is to have the generating 
and sensing coils move with respect to each other, such as is taught in United States Patent No. 3,644,825, 
issued February 22 1972 to Paul D. Davis, Jr., entitled Magnetic Detection System for Detecting Movement 
of an Object Utilizing Signals Derived from Two Orthogonal Pickup Coils. The motion of the coils generates 
changes in the magnetic field, and the resulting signals then may be used to determine the direction of the 

20 movement or the relative position of the generating and sensing coils. While such an approach removes 
some ambiguity about the position on the basis of the field sensed, its accuracy is dependent on the 
relative motion, and it cannot be used at all without the relative motion. 

Another approach that has been suggested to provide the additional required information is to make 
the magnetic field rotate as taught in Kalmus, "A New Guiding and Tracking System," IRE Transactions on 

25 Aerospace and Navigational Electronics, March 1962, pages 7 to 10. To determine the distance between a 
generating and a sensing coil accurately, this approach requires that the relative orientation of the coils be 
maintained constant. It therefore cannot be used to determine both the relative translation and relative 
orientation of the generating and sensing coils. 

United States Patent No. 3,868,565, issued February 25, 1975, to Jack Kuipers entitled Object Tracking 

30 and Orientation Determination Means, System and Process teaches a tracking system for continuously 
determining at the origin of a reference coordinate system the relative translation and orientation of a 
remote object. The tracking system includes radiating and sensing antenna arrays each having three 
orthogonally positioned loops. Properly controlled excitation of the radiating antenna array allows the 
instantaneous composite radiated electromagnetic field to be equivalent to that of a single loop antenna 

35 oriented in any desired direction. Further control of the excitation causes the radiated field to nutate about 
an axis denoted as a pointing vector. This tracking system is operated as a closed-loop system with a 
computer controlling the radiated-field orientation and interpreting the measurements made at the sensing 
antenna array. That is, an information feedback loop from the sensing antenna array to the radiating 
antenna array provides information for pointing the nutation axis toward the sensing antenna array. 

40 Accordingly, the pointing vector gives the direction to the sensing antenna array from the radiating 
antenna array. The proper orientation of the pointing vector is necessary for computation of the orientation 
of the remote object The signals detected at the sensing antenna include a nutation component. The 
nutating field produces a different nutation component in each of the three detected signals. The 
orientation of the sensing antenna array relative to the radiated signals is determined from the magnitudes 

45 of these components. 

United States Patent No. 4,054,881, issued October 18, 1977 to Frederick H. Raab, entitled Remote 
Object Position Locater, describes apparatus for determining the position and orientation of a remote 
object relative to a reference coordinate frame comprising: a plurality of radiating means having 
orthogonal components centered about the origin of the reference coordinate frame; means for applying to 

so the plurality of radiating means electrical signals which generate a plurality of electromagnetic fields, the 
electromagnetic fields being distinguishable from one another; a plurality of receiving means disposed on 
said remote object, said receiving means having orthogonal components for detecting and measuring 
components of the electromagnetic fields; and analyzing means associated with the receiving means for 
converting the received components of the electromagnetic fields into remote object position and 

55 orientation relative to the reference coordinate frame with at most two ambiguous combinations of 
orientation and position, the analyzing means operating open loop with respect to the radiating means. 
This specification thus teaches a magnetic or near-field non-tracking system for determining, at a remote 
object, the position of the remote object with respect to a reference coordinate system. The orientation of 
the remote object can be determined, at the remote object, with respect to the reference coordinate system 

so by using an iterative computational scheme. This is accomplished by applying electrical signals to each of 
three mutually orthogonal radiating antennas, the electrical signals being multiplexed with respect to each 
other and containing information characterizing the polarity and magnetic moment of the radiated 
electromagnetic fields. The radiated fields are detected and measured by three mutually orthogonal 
receiving antennas, having a known spatial relationship to the remote object, which produces nine 

65 parameters. These nine parameters, in combination with one known position or orientation parameter are 
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sufficient to determine the position and orientation parameters of the receiving antennas with respect to 
the position and orientation of the radiating antennas. 

It has now been recognised that remote object position and orientation can be calculated with a 2 — 3 or 
3 — 2 arrangement rather than the 3 — 3 arrangement described in the earlier specification. Thus, according 
to the present invention, apparatus for determining the position and orientation of a remote object relative 
to a reference coordinate frame and comprising the elements set out above is characterized by the features 
listed in the characterizing clause of Claim 1. 

Thus, by means of the invention, either the transmitter or the receiver is reduced from three 
components to two components and this significantly reduces the cost and complexity of the apparatus. 

The invention is thus directed to a technique for determining relative position and orientation between 
an electromagnetic source and sensor which employs far-field coupling between the source and sensor. In 
the case of large separation distances between the source and sensor, the present invention presents the 
advantages of smaller size, weight and power requirements with longer range and minimization of 
multipath signals (field distortion) by using short pulse transmission. Transmission and processing speed 
1S is increased by applying electrical signals to the radiating means that consist of only two discrete excitation 
states. The transmitted electromagnetic fields are multiplexed and thus are distinguishable from one 
another. 

Processing speed in the present invention is substantially increased over that of the prior art and the 
cost and complexity of either the transmitter or receiver is substantially reduced by providing one of the 
plurality of radiating means and the plurality of receiving means with only two orthogonal components and 
the other with three orthogonal components. The analyzing means associated with the receiving means 
converts the received components of the transmitted electromagnetic field into remote object position and 
orientation relative to the reference coordinate frame with at most two ambiguous combinations of 
orientation and position. 

25 Further advantages of the invention and of the features of the sub-claims will become apparent from 
the following description of one form of apparatus embodying the invention, namely a landing aid system, 
and of a number of possible modifications or alternatives which may be incorporated in the system. The 
system which constitutes an example only of the invention will be described with reference to the 
accompanying drawings, in which: 

30 Figure 1 is a partly block, side elevationai view of a landing aid system not forming part of the invention 
in which three components are provided both on the transmitter and receiver; 

Figure 2 is a graphical representation showing the relationship between electric field strength 
(ordinate) and distance from a radiator (abscissa); 

Figure 3 is a simplified representation of an electric field associated with a current-carrying electric 

3$ dipole; 

Figure 4 is a graphical representation of the location coordinate system of the remote object with 
respect to the location of the origin of the reference coordinate frame; 

Figure 5 is a graphical representation of the orientation coordinate system of the remote object with 
respect to the reference coordinate frame; 
40 Figure 6 is a graphical representation of the amplitude of the signals applied to the transmitting 
antennas, with respect to time, in the case where the signals are frequency division multiplexed; 

Figure 7 is a block diagram of a portion of the receiver of the system; 

Figure 8 is a graphical representation of the far-field electromagnetic coupling of a three axis sensor to 
a three axis source; 

45 Figure 9 is a graphical representation of the far-field electromagnetic coupling of a three axis sensor 
with a three axis source of unknown orientation and structure; 

Figure 10 is a flow chart for the computations carried out in a three-state power solution for remote 
object position and orientation; 

Figure 11 is a flow chart of the computations carried out in a two-state power and dot product solution 
so for remote object position and orientation; 

Figure 12 is a graphical representation of the signals applied to the transmitting antennas, with respect 
to time, in the case where the signals are time division multiplexed; 

Figure 13 is a schematic representation of a transmitter employed in a time division multiplexed 
system; and 

55 Figure 14 is a schematic representation of a transmitter employed in a frequency division multiplexed 
system. 



Description of the preferred embodiments 
Apparatus 

eo Three axis transmission and three axis sensing with frequency division multiplexing 

Although the invention may have utility in a number of environments, a long distance landing system 
is described in detail. Referring to Figure 1, which does not form part of the invention but is described first 
to aid an understanding of the invention, a landing aid system 10 includes ground based components 30 for 
radiating an electromagnetic field and airborne components 20 for receiving the electromagnetic field and 

•55 determining the position and orientation of the airborne components 20 with respect to the ground based 
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components 30. The ground based components include a signal generator 31 coupled in parallel to power 
amplifiers 32, 33 and 34. A ground antenna array 40 includes orthogonal electric dipole antennas 41, 42, 
and 43 (denoted X f Y f Z) coupled to power amplifiers 32, 33, and 34, respectively. The dipole antennas 41, 
42 and 43 are short relative to the wave length of the carrier frequency so that they each produce an electric 

5 dipoie-field pattern unique to each antenna. A monitor receiver 44 is coupled to signal generator 31, spaced 
from ground antenna array 40 and has an orthogonal antenna array 45 for receiving electromagnetic 
radiation from ground antenna array 40. The separation distance of monitor receiver 44 from the ground 
antenna array 40 is such that the electromagnetic field has a far-field component substantially in excess of 
the near-field component. Monitor receiver 44 provides a means of verifying the electromagnetic 

to transmission from ground antenna array 40. Airborne components 20 include the sequential coupling of an 
orthogonal receiving antenna array 21, a signal amplifier group 52, a frequency translator group 53, a 
signal processor group 54, a computer 50 and a display 51. More specifically, the antenna array 21 includes 
receiving dipole antennas 22, 23 and 24 (denoted U, V, W) coupled sequentially to signal amplifiers 25, 26 
and 27, respectively, frequency translators 55, 56 and 57, respectively, and signal processors 58, 59 and 60 

75 respectively. 

The landing aid system 10 operates "open loop" in that the only communication between the airborne 
components 20 and the ground based components 30 is the radiated electromagnetic field from the ground 
based components 30. There need be no communication from the airborne components 20 to the ground 
based components 30 in order to establish the position and orientation of the receiving antenna array 21 

20 with respect to the ground antenna array 40. Further, the landing aid system 10 allows simultaneous use by 
any number of remote users. In addition to providing the capability for measuring position and orientation, 
the signals radiated by the ground antenna array 40 can provide a one-way data link from the ground based 
components 30 to the receiving antenna array 21. The link can carry information such as transmitter 
identification, transmitter power, field distortion corrections, locations of nearby obstacles, the location of 

25 the landing site relative to the ground antenna array 40 and wind direction. 

Referring to Figure 2, the field produced by excitation of a dipole antenna can be separated into two 
components referred to as the near-field and the far-field components. In use of the present invention, the 
separation distance of the remote object from the transmitting means is limited to far-field conditions. The 
far-field component of the transmitted electromagnetic radiation decreases as the inverse of distance 

30 between the remote object and the transmitter increases. The intensity of the far-field depends on the 
relative size of the antenna and the wave length of the excitation frequency. For electrically short antennas, 
as the wave length of the excitation frequency is shortened, or the excitation frequency is increased, the 
strength of the far-field component increases. The far-field component of electromagnetic radiation is 
generally used for long distance communications and navigation. On the other hand, the near-field 

3S component of electromagnetic radiation decreases with the cube of the distance from the antenna 
preventing its detection at large distances. The intensity of the near-field is not a function of frequency and 
it can be quite high at short distances. When using the far field component some additional field distortion 
occurs because of surrounding objects. The amount of distortion resulting from surrounding objects 
depends on the conductivity and permeability of these objects and their size and location relative to the 

40 receiving and transmitting antennas. It is possible to predict and compensate the distortion caused by 
nearby fixed objects and hence essentially remove position and orientation errors caused by these objects. 

The ground based components 30 generate a far-field landing aid signal. The signal generator 31 
generates the electrical signals which excite each of the antennas 41, 42 and 43. The signal must be 
multiplexed so that the receiving antenna array 21 can distinguish the electromagnetic radiation from each 

45 of the antennas 41, 42 and 43. Although the list is not exhaustive, the electromagnetic radiation transmitted 
from each of the antennas 41, 42 and 43 may be distinguished by using time division multiplexing, 
frequency multiplexing, phase multiplexing and spread spectrum multiplexing. Additionally, the electrical 
signal may contain information characterizing the phase of the electromagnetic radiation. A simple 
example would be to include a timing pulse whenever the signal goes positive. Alternatively, if frequency 

so multiplexing is used, the excitation to each of the antennas 41, 42 and 43 is advantageously coherent That 
is, periodically all of the signals go positive simultaneously (see Figure 6). Additionally, the data frequency 
determines the spacing between the carrier frequencies, and is thus the basic reference frequency of signal 
generator 31 . The data frequency is labeled f 0 in Figure 6. Advantageously, the reference frequency will be 
derived from a temperature compensated crystal oscillator in the 10 MHz range and frequency selection 

55 will be in 10 kHz steps. 

The three power amplifiers 32, 33 and 34 boost the outputs of the signal generator 31 to a level 
sufficient to produce the desired power with the given antenna. To make efficient use of the power 
available, a switching power may be used. For example, either class D (carrier frequency switching) or a 
class S (high frequency switching) modulator can be used. An RFI filter is advantageously also included. 

60 The ground antenna array 40 includes mutually orthogonal dipole antennas 41, 42 and 43 and may be 
located near the landing pad. The relationship of the landing pad to the ground antenna array 40 can be 
included in the one-way data stream to the airborne components 20. Antenna design in the far-field context 
is highly dependent on the operating carrier frequency. For a long distance landing aid system 10, a carrier 
frequency of 220 MHz is appropriate. A dipole antenna whose length is approximately one-tenth of the 

55 wave length of the carrier frequency would give a dipole length of approximately 12 centimetres. 
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The monitor receiver 44 is similar to an airborne receiver, but omits position/orientation computations, 
data decoding, and display. Its function is to ensure that electromagnetic field amplitudes and phases 
radiated from the ground antenna array 40 are correct. When deviations are found, change instructions are 
issued to the signal generator. If signals cannot be maintained within prescribed tolerances, the monitor 

5 can place an out-of-tolerance message in the data stream. Of course, it can be appreciated that the monitor 
receiver 44 is not necessary to an embodiment of this invention. 

The airborne components 20 of the landing aid system 10 for a frequency division multiplexing 
embodiment are shown in Figure 1 and provide a separate signal processing path for each of the signals 
from the receiving dipole antennas 22, 23 and 24. Although there is additional discussion below of various 

70 multiplexing alternatives, it can readily be appreciated that if time division multiplexing were used, a single 
path could be switched among antennas 22, 23 and 24. 

Figure 7 shows a more detailed block diagram of a signal path, in particular the U antenna signal path. 
For practical reasons, the amplifier group 52 and the frequency translator group 53 are broken into several 
components and spread throughout the signal path. More specifically, the U amplifier 25 of the amplifier 

is group 52 includes a preamplifier 61, gain control 62, amplifier 63 and amplifier/clipper 64. The frequency 
translator 55 of the frequency translator group 53 includes a band pass filter 65, mixer 66, second band pass 
filter 65', mixer 67, low pass filter 68 and synthesizer 69. Accordingly, the receiving dipole antenna 22 is 
coupled to the preamplifier 61, band pass filter 65, gain control 62, amplifier/clipper 63, mixer 66, second 
band pass filter 65', amplifier/clipper 64, mixer 67 and low pass filter 68. Synthesizer 69 is connected to the 

20 signal processor 58. 

The signal processor 58 includes parallel combinations of the sequential connection of a mixer, an 
integrator and a sample-and-hold block. More specifically, each of the parallel paths has a mixer 70 to 75, 
an integrator 76 to 81 and sample-and-hold blocks 82 to 87. The outputs from the sample-and-hold blocks 
82 to 87 are coupled to the computer 50 and in turn coupled to the display 51. In this embodiment, there are 

25 six parallel paths for processing the signal from the U antenna. There is one path for signals received from 
the Y transmitting antenna received by the U receiving antenna, denoted Y/U. Similarly, there is one path 
for signals received from the X transmitting antenna received by the U receiving antenna, denoted X/U. The 
transmitted signal from the Z antenna has two frequencies for carrying a binary code and requires two 
paths, denoted Z,/U and Zq/U. Further, during acquisition two additional signals are processed so data 

30 transmitted by the Z antenna is not lost. These signal paths are denoted Z^Uq and Z^Uq and have 
negligible output when the receiver is locked on the transmitted frequencies. 

The metal aircraft upon which the receiving antenna array 21 is mounted causes some distortion of the 
electromagnetic fields received by the antenna. Unless the aircraft is very close to the transmitter, this 
distortion may be described by a linear transformation which maps the free space fields into the three 

35 antennas 22, 23 and 24. For example, a field aligned exactly with the length of the aircraft will also appear in 
the transverse and vertical receiving antennas. This effect is constant for a given aircraft and installation. It 
is easily corrected by applying an inverse linear transformation to the measured data. 

The input bandwidth of the amplifier group 52 is advantageously restricted to the 219—221 MHz band 
after which the signals are boosted to a suitable level and noise impulses are clipped. Accurate gain control 

40 is used to obtain maximum effectiveness in clipping noise. Alternatively, an impulse detector may be 
applied at this point to shut off amplifier 25 when an impulse occurs. The signal is now translated 
downward to a convenient intermediate frequency such as 10 MHz. The bandwidth is also reduced to 100 
kHz. After the final amplification and clipping, the signal is translated downward to approximately 100 kHz 
for final processing. 

45 The mixing frequencies required to accomplish the necessary frequency translation are synthesized by 
standard techniques. The first mixing frequency is selectable in 10 kHz steps from 209—211 MHz. This 
allows any selected signal in the 219—221 MHz band to be translated to 10 MHz. The second mixing 
frequency of 10.1 MHz is fixed and translates the 10 MHz intermediate frequency to the 100 kHz processing 
frequency. For initial acquisition, these frequencies are synthesized from a stable reference oscillator. After 

so acquisition, they can be locked to the received signal to remove any frequency error. 

Signal processor group 54 must acquire the received signals, establish the timing reference, make 
measurements for position/orientation computations, and decode transmitted data. To do this, it uses a 
series of phase-locked loops, frequency dividers, and integrators. Interface with computer 50 is 
accomplished by an A/D converter and a suitable buffer. Signal acquisition is accomplished by the 

55 equivalent of a pair of phase-locked loops. In this particular example, as stated, frequency multiplexing is 
used and data is carried only by one radiating antenna (the Z signal shown in Figure 6). The frequencies for 
carrying data, i.e. ones and zeros, on the Z signal are referred to as the mark and space frequencies. 
Accordingly, the phase-locked loops can operate at, for example 110 KHz to correspond to the mark and 
space frequencies of the Z signal. Loop bandwidth may be changed for initial acquisitions and later 

oo tracking, but in either case, it will be sufficiently low to cause the loop to ignore the effects of the frequency 
shift keying. The 10 kHz reference timing is obtained as the difference between the frequencies of the two 
oscillators generating the mark and space frequencies. Actual implementation can use measurements of 
both the sine and cosine integrals for measurement and locking, respectively. 

Signal measurements are made by mixing a received signal with a locally generated signal and 

55 integrating the product. A coherent set of mixing frequencies (for example, 120, 110, 100, 90 kHz) 
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corresponding to the four transmitting frequencies is synthesized from the 10 kHz reference frequency. The 
integrators 76 to 81 are advantageously reset about every .001 seconds by the reference signal. The value 
in each of the integrators is transferred to the sample-and-hold circuits 82 to 87 just prior to the resetting of 
the integrators 76 to 81. 

5 Decoding of the data and averaging of the measurements is accomplished by software. The computer 
50 can measure signal amplitude and signal-to-noise ratio on a sample-by-sample basis. Navigation 
measurements of the X and Y signals are accomplished simply by summing an appropriate number of .001 
second samples. A similar procedure is used on the Z channel for initial acquisition. When measurements 
indicate a satisfactory signal-to-noise ratio, data may be extracted by comparing the Z-mark samples to the 

to Z-space samples. Z-navigation information is based on an average of those samples corresponding to the 
data received. That is, only the Z-mark or Z-space sample at a given sampling point is used, depending on 
the decision about which carrier was transmitted during that interval. 

The computer and display can be common to both the long distance landing aid system 10 and the 
landing aid system disclosed in United States Patent No. 4,054,881 to Frederick H. Raab issued October 18, 

is 1977 entitled Remote Object Position Locater. This is particularly advantageous for reducing cost and for 
simplification of equipment. Further, an aircraft may use the present long distance landing aid system 10 to 
navigate to within a few kilometres of the landing point and then acquire signals from the landing aid 
systems disclosed in the aforementioned patent for final approach guidance. The computer and display can 
be anything suitable and are therefore not discussed here in detail. 

20 

Two axis transmission or two axis sensing with frequency division multiplexing: 

Although Figures 1 and 7 detail a landing aid system 10 utilizing three transmitting antennas, 41, 42 
and 43, and three receiving antennas 22, 23 and 24, in accordance with an embodiment of the invention a 
landing aid system utilizing two transmitting antennas 41 and 42, and three receiving antennas 22, 23 and 

25 24, or a landing aid system utilizing three transmitting antennas 41 , 42 and 42, and two receiving antennas 
22 and 23, may be provided. Two axis transmission with three axis sensing simplifies the transmitter. This 
arrangement also brings about an increase in processing in the case where time division multiplexing is 
used to distinguish the signals applied to each axis of the transmitting antenna array. Three axis 
transmission with two axis sensing simplifies the receiver. However, the use of two antennas for either 

30 transmitting or receiving does add an additional ambiguity to the system. This can be corrected by, in 
addition to specifying that the airplane is flying right-side-up or upside down, specifying that the airplane is 
approaching the landing site from either the north or south, or the east or west. Apparatus for transmitting 
with only two orthogonal antennas is the same as that previously described with reference to Figure 1 
except that only two of the X, Y or Z signals paths are necessary. Apparatus for receiving these signals is 

35 the same as that previously described with reference to Figures 1 and 7 except that the remaining signal 
paths contain fewer parallel paths since the signal from one of the X, Y or Z transmitting antennas is not 
present. Apparatus for receiving three transmitted signals with only two receiving antennas is the same as 
that previously described with respect to Figures 1 and 7 except that only two signal paths are necessary for 
the two orthogonal receiving antennas. 

40 

Time division multiplexing: 

Figure 12 depicts a pulsed carrier wave signal format suitable for use in a time division multiplexed 

system. The three axes of the transmitting antenna are excited sequentially by signals of the same 

frequency. The durations of the three pulses are known (fixed), with the X-axis excitation pulse longer than 
45 the others to allow the receiver to establish synchronization, thereby knowing which received signals to 

attribute to which transmitting axis. 

To allow rejection of multipath effects, a "dead space" might be inserted between the pulses to allow 

time for echos to die out. If multipath interference were no problem, all three axes could be excitated 

simultaneously by signals of different frequencies or by signals modulated by different spread spectrum 
so codes. These are ordinary engineering design decisions that must be made for each application of the 

disclosed concept. 

Formats for two-state excitation are similar but simply omit excitation of one axis. If two-axis reception 
is used, a three-state excitation pattern as described above is still required. 

The carrier frequency for these signals would normally be in the range of 300 to 3000 MHz with present 
55 technology. The excitation pattern could be repeated at frequencies in the range of 1 kHz to 30 MHz. 
Rgure 13 depicts a block diagram of a transmitter for a time division multiplexed system. Note that for 
two-state transmission, the Z-axis of the antenna and the associated driving circuitry is omitted. 

All signals in the transmitter are derived from a stable oscillator 200 by a frequency synthesizer 201. 
The derived radio frequency signals are switched to the power amplifiers 207, 208 and 209 by gates 203, 
60 204, and 205, which operate under the control of a sequencer 206. The power amplifiers 207, 208 and 209 
produce excitation voltages w Xf w y , and w z as inputs to the antenna axes 210, 211, and 212, respectively. 
The antennas 210, 21 1 and 212 are dipoles that are short relative to the wavelength of the carrier frequency 
so as to produce an electric dipole-field pattern. 

A receiver suitable for use in a time division multiplexed system is illustrated in Figure 14. Signals are 
55 received by short dipole antennas 213, 214, and 215 and preamplified by preamplifiers 216, 217, and 218. 
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For two-axis reception of three-state transmissions, one of the receiving antennas and the circuitry 
associated therewith can be omitted. 

After preamplrfication, the three received signals are converted to an intermediate frequency by mixers 
219, 220 and 221, which are driven by signals produced by an oscillator 222, and a synthesizer .223. Note 
that all signals and timing in the receiver are derived from one master oscillator. Not shown is apparatus for 
phase-locking to the received signal, which may be added and is standard technology. 

Intermediate-frequency signals are amplified by amplifiers 224, 225, and 226. The amplified 
intermediate-frequency signals are mixed with signals of same frequency in mixers 227, 228, and 229. The 
outputs of these mixers are integrated by integrators 230, 231, and 232, and sampled by sample-and-hold 
circuits 233, 234, and 235; outputs are acquired by the computer 239, which performs the required 
mathematical operations to extract position and orientation information, which is displayed by a display 
240. The mathematical operations hereinafter developed are equally applicable to time division 
multiplexed and frequency division multiplexed signal formats as well as many other types of signal 
formats. 

Operation 

Referring now back to the frequency division multiplexing embodiment illustrated in Figures 1, 6 and 7, 
if unambiguous measurements are desired geometrical considerations result in inclusion of a timing 
reference in the transmitted signal. They also require airborne components 20 to measure the signal 
components induced in each receiving antenna 22, 23 and 24 by each transmitting antenna 41, 42 and 43. 
These requirements and any additionally desired data transmission form the constraints on signal format. 
While any choices are possible, coherent frequency division multiplex/frequency shift keying may be 
suitable for many general purpose users. 

It should be noted that in order to facilitate the orderly development of a position and orientation 
finding algorithm, the three transmitting antennas will be designated a three axis source 98, and the three 
receiving antennas will be designated a three axis sensor 100. 

Far field coupling: 

Excitation of an electric dipole or loop (magnetic dipole) antenna produces terms that vary as 1/p , 1/p , 
30 and 1/p, which are referred to as quasi-static (near) field, induction field, and far field, respectively. In the 
invention a far field coupling is used. At large distances (p»A/2n), the far-field terms dominate and the 
resultant electric and magnetic fields form essentially plane waves. The electric and magnetic field vectors 
are orthogonal to each other and both are othogonal to the direction of propagation. The cross product of 
the electric and magnetic field vectors, called the Poynting vector, represents power flow, and is oriented in 
35 the direction of propagation. 

It is convenient to think of far-field coupling in terms of the behaviour of electric dipoles, although 
essentially the same relationships hold for magnetic dipole (loop) sources and sensors. The electric field 
resulting from the excitation of an electrically short dipole is: 
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where the excitation current is I cos cot, the antenna length is I, and A is the wavelength of the carrier 
45 frequency. The off-axis angle 6 and the field pattern providing for each antenna, as defined by Equation (1 ) 
are shown in Figure 3. Note that in contrast to the near-field, the far-field intensity varies as the inverse of 
distance and is frequency dependent. 

The magnetic field vector is related to the electric field vector by the free-space impedance n— 377Q, 

thus 



|E| 

| H |- < 2 > 



In the receiving or (or sensing) mode, a dipole has the same pattern as it does in transmitting. 

An elementary dipole sensor therefore produces an output proportional to the cosine of the angle 
between the electric field vector and the dipole axis. Note that variation of the field strength with the sine of 
the off-axis angle 5 is a characteristic of elementary short dipoles. This simple variation does not apply to 
60 dipoles whose lengths are a significant portion of a wavelength (0.1\ or greater) or arrays of coupled 
colinear elements. For example, the field strength produced by a half-wave dipole varies as 

n 

cos ( — cos 5)/sin 8. 
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A coordinate system for determination of the position of the receiver relative to the transmitter is 
shown in Figure 4. The X, Y and Z axes are aligned with north, east, and vertical, respectively, and are 
centered at the centre of transmitting ground antenna array 40. Location of the receiving antenna array 21 
may be specified in either rectangular coordinates (x, y, z) or in polar coordinates (a, P, p). It may also be 
5 specified by the distance p and two of the three orientation angles i(jx, ipy, or ipz. 

Measurement of the three transmitted signals from ground antenna array 40 as received by the set of 
three orthogonal receiving antennas 22, 23 and 24 produces nine parameters which are sufficient to 
determine the six position and orientation parameters. As noted earlier this assumes one orientation or 
position parameter is independently determined. While there are a variety of computational algorithms 
to that can be used, it is conceptually easiest to begin by using relative amplitudes to determine position. 

The first step in synthesizing a position and orientation finding algorithm is the definition of 
coordinates and vector-matrix formulations relating sensor output to source excitation. 

The geometric relationship between the three-axis source 98 and the three-axis sensor 100 is shown in 
Figure 5. The source coordinate frame X t -V, -Z, is defined by the axes of the source 98. Alignment of 
ts these axes with some convenient natural reference such as north, east and down is assumed. The source 
axes can be effectively aligned with any desired coordinate frame by altering the excitation. Similarly, 
coordinates measured in the source coordinate frame can be converted to any desired coordinate frame 
mathematically. 

The sensor position is specified in rectangular (x, y, z) or spherical (a, 0, p) coordinates defined relative 

20 to the source coordinate frame. Sensor orientation is specified by a sequence of three rotations. Azimuthal 
rotation by first turns the sensor about its Z axis from +X toward +Y. The elevation rotation by 6 then 
turns the sensor about its Y axis from +X to -Z. Finally, a roll rotation by <p turns the sensor about its X axis 
from +Y to +Z. Note that in the zero-orientation condition, the three sensor axes are parallel to the 
corresponding source axes, and that the order of the rotations cannot be interchanged without changing 

2s the values of ip, 6, and (p. 

The excitation of a three-axis electric dipole source 98 and the resultant three-axis sensor output are 
most conveniently described in vector notation. The excitation of the source is therefore represented by 
i=[fx/ fiy* fiJ T - The lengths of the three dipoles are assumed to be identical, hence f 1x , f 1y , and f lz represent 
the amplitudes of the currents exciting the dipoles of X-axis, Y-axis, and Z-axis orientation, respectively. 

30 Now let the output of a three-axis sensor be similarly represented by^=[f 3x , f 3y , faJ T , and consider the 
coupling between that sensor and a similarly aligned source £. Figure 8 depicts a three-axis source 102 and 
a three-axis sensor 1 00 whose coordinate frames are aligned. Since the sensor 1 00 is located on the X 2 axis, 
the sensor 100 is in the null of the X 2 dipole, hence the X 2 excitation produces no sensor response on any 
axis. The source Y 2 axis is parallel to the sensor Y 3 axis, and therefore produces a response in that axis. 

35 However, the electric field resulting from Y 2 axis excitation is orthogonal to yye sensor Z 3 axis and hence 
produces no Z 3 response. Coupling between the source Z 2 excitation and the sensor Y 3 and Z 3 axes is 
analogous. 

If the three-axis source excitation is represented as a vector £ and the three-axis sensor output is 
similarly represented as a vector the source-to-sensor coupling can be described by: 



40 



45 



55 
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i3= — s £=- 



0 0 0 
0 1 0 
0 0 1 



(3) 



The factor C accounts for excitation and sensing constants common to all axes. Note that the far-field 
coupling matrix S (shown above) is degenerate and differs from the near-field coupling matrix. 

The coupling between a source 98 and sensor 100 of arbitrary position and orientation (Figure 5) can be 
so determined by inserting orthogonal rotation matrices into Equation (3). These matrices are based upon 
position azimuth and elevation (a and P) and orientation azimuth, elevation, and roll (ip, 9, and as shown 
in Table 1. Note that the subscript defines both the type of transformation and its independent variable. 
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TABLE 1 



Position 



Orientation 



Azimuth 


~ cos a sin a 


Q 




"~ rmc ill 


Qin ill 

Oil 1 


0 1 


rotates 














X into Y T = 


-sin a cos a 


0 


I»= 


—sin tp 


cos 


0 




0 0 


1 




0 


0 


1 


Elevation 


" cos P 0 


-sin p 




cos 6 


0 


—sin 0 


rotates 














X into -Z Tp= 


0 1 


0 




0 


1 


0 




sin p 0 


cos p 




_sin 0 


0 


cos 0 


Roll 








1 


0 


0 


rotates 












sin $ 


Y into Z 


not applicable 


I»- 


0 


cos $ 










0 


-sin <{> 


cos <p 



Inverses 



ri-X-- Op D-^I-o I- 



— Consider first the coupling between the source and a zero-orientation sensor (whose output isi*), 
located at (a P, p), as shown in Figure 5. The excitation fa of an equivalent source 102 whose X-axis is 
aligned with the line connecting the source 98 and sensor 100 can be determined by rotating the excitation 
vector of the real source 98 by position azimuth and elevation, thus 



(4) 



The coupling to a similarly aligned equivalent sensor £ then has the same form as Equation (3), i.e., 
i3=(C/p) §£. The output of the zero-orientation sensor is then found by applying inverse position rotations, 
thus 



C C 

&= I-a I- P & Ip la i= Q £ 

P P 

The equivalent sources and sensors used above are listed in Table 2. 

TABLE 2 



(5) 



50 



55 



Symbol 


Name 


Definition 


fi 


" Source 




f 2 


Position-frame source 




f 3 


Position-frame sensor 


c 
p 


u 


Zero-orientation sensor 






Sensor 


is = Xt> Ib I* i» 
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Utilizing Table 1 and equations (3) and (4) f, can be expanded as follows: 



A= §lplai 

P 



(6) 



w 



75 



0 0 0 "ITcos (3 0 -sin p 
0 10 0 1 0 
_ 0 0 1 JLsin P 0 cos P J _ 
0 0 0 

—sin a cos a 0 

_ cos d sin p sin a sin (3 cos P _ 



cos a sin a 0 
-sin a cos a 0 
0 0 1 



(7) 



(8) 



The row of zeros in Equation (8) implies that no source excitation can produce a radial (position-frame 
20 X 3 ) component. 

A fixed, three-state excitation pattern based upon the source axes is given by: 



(9) 



This is the same excitation pattern used by the near-field large-angle algorithm. The fields produced at 
30 the sensor location in response to these excitation vectors are then: 





" r 




"0" 




" 0 * 


25 A (S1)= 


0 

-'0. 


A <S2H 


1 

.0. 


A (S3>- 


0 

. 1 _ 



A (S1=- 



35 



40 



A (S2)=- 



0 

-sin a 
_ cos a sin p 
0 

cos a 

sin a sin p 



C 
P 



-y/Vx^+y* 

_x z/pVx 2 +y 2 , 
0 

x/Vx 5 ^ 
-y z/pVx^+y 2 _ 



(10) 



(11) 



and 



45 



A (S3)=- 



50 





" 0 




0 


c 




C 






0 




0 


p 




p 


^-vV+y^p 




_cos p 





(12) 



The output of a three-axis sensor of arbitrary orientation (ip, 9, (p) is determined by applying orientation 
azimuth, elevation, and roll rotations to the output of the equivalent zero-orientation sensor, thus: 



55 



A— T A Sa Xp Ja Xn I- a I-p § Ip la A 

P 



(13) 
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Position determination: 

Since the sensor orientation is unknown at this point in the processing, use must be made of 
orientation-independent signal parameters. Three such parameters are: 

1. Signal power, obtained by dot products by sensor response vectors with themselves; 

2. Dot products between different sensor response vectors, analogous to the angles between the response 
vectors; and 
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3. The amplitudes of the cross products between different sensor responses, analogous to the 
perpendiculars to the planes defined by those responses. 

Several algorithms are derived subsequently; the choice of algorithm depends upon the application. 

5 

1. Three-state power solution: 

The sensor responses from all three states of source excitation Equations (10), (11), and (12) can be 
converted to received power and processed to yield position in a manner similar to that used for near-field 
applications. The three "power" responses are obtained by taking dot products of the three sensor 
to response vectors with themselves. Since the sensor orientation is determined by a set of orthogonal 
rotations, the power is invariant under sensor rotation. The three power outputs are: 

C 2 

P(Si(= [sin 2 a+cos 2 a sin 2 P] (14) 

c 2 

P(S2|= [cos 2 a+sin 2 a sin 2 0] (15) 

P 2 



20 



45 



P(S3>= cos 2 p (16) 



25 Distance p is obtained from the sum of the three powers, which is independent of a and P, thus the 
calculated distance p equals: 



f 2C 2 

30 \ P(Sl)-t-P(S2)+P(S3) 



Rearrangement of Equation (16) then yields: 



4 



P 2 



35 I P |=arccos v I P(S3) (18) 

C 2 

Substitution of p and | p| into equation (15) then produces: 



40 I p 2 



P(S2)-1 



ct A =arcsm (19) 

sin 2 | p| -1 

The position defined by a A , | P | , and p contains an eight quadrant ambiguity (as in the near-field 
algorithm), which is reduced to a two-quadrant ambiguity by the signs of the dot products. The three 
possible dot-products are: 

so C 2 

v(S1, S2)= (-sin a cos a+sin a cos a sin 2 P) 

P 2 

C 2 

55 = (sin 2a) (cos 2 p) (20) 

2p 2 

C 2 

v(S2, S3)= 

60 p 2 

C 2 

(sin a sin P cos P)= sin a sin 2P (21) 

2p 2 

65 and 
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v(S3, S1)=- 



(cos a sin P cos (3)=- 



C 2 



2P S 



cos a sin 2(3 



(22) 



Inspection of Table 3 shows that the polarities of any two of these dot products reduces the quadrant 
to ambiguity from 8 to 2. 

TABLE 3 











Position coordinates 




Dot products 




Quadrant 


X 


y 


z 


a 


P 


v(S1, S2) 


v(S2, S3) 


v(S3, SI) 


1 


+ 


+ 




0°... 90° 


0°...-h90° 




+ 


+ 


2 




+ 




90°... 180° 


0°...+90° 


+ 






3 








-90°... -180° 


0°...+90° 








_ 4 


+ 






0°... -90° 


0°...+90° 


+ 




+ 


5 


+ 


+ 


+ 


0°... 90° 


0°...-90° 








6 




+ 


+ 


90°... 180° 


0°...-90° 


+ 




+ 


7 








-90°...-180° 


0°...-90° 




+ 


+ 


8 


+ 






0°...-90° 


0°...-90° 
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A somewhat more direct solution can be obtained if the power responses are formulated in rectangular 
coordinates: 
First, 



P(S3) 



x 2 +y 2 



C 2 



(p 2 -z*) 



Geometrical Similarity then requires that 

C 2 y^z 2 



(23) 



P(S1)= 



(p'-x 2 ) 



(24) 



50 



and 



P(S2)= 



xHz 2 



C 2 

^4 



(p'-y 2 ) 



(25) 



55 Range p is first found by using Equation (17). Values of x 2 , y 2 , and z 2 are then found by substitution of 
measured "power" and p 2 into Equations (22), (24), and (25). Figure 10 illustrates a flow diagram for the 
computations involved in the implementation of a three-state power solution for position. Figure 10 also 
goes on to illustrate the flow diagram for the computations involved in calculating orientation from the 
three-state power solution for position. The mathematical operations for calculating orientation are 

eo presented later. 

2. Two-state power and dot product solution: 

A two-state large-angle algorithm similar to that for near-field operation can be developed by using the 
normalization 

65 
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Cx Cy Cz 

X= Y= Z= (26) 

P 2 P 2 P 2 

Equations (24) and (25) then become: 

P(S1)=Y 2 +Z 2 (27) 

and 

P(S2)=X 2 +Z 2 (28) 
The first dot-product Equation (20) is then also converted to normalized rectangular coordinates 

C 2 , x w y v # VS^+y 5 " % 2 

= -X Y (30) 



v(S1 r S2)= (___)(___)( ) (29) 



" If the dot-product Equation (30) is zero, either X=0 or Y=0 or both; which of these is the case is readily 
determined by whether P(S2) >P(S1), P(S1) >P(S2) or P(S1)=P(S2), respectively. If v(S1, S2)*0, then 
equation (30) can be rearranged into: 

25 -v(S1, S2) 

Y = (31) 

X 



The difference of (28) and (27) then eliminates Z,: 

P(S2)-P(S1)=X 2 -Y 2 (32) 
and substitution of Equation (31) produces an equation containing X 2 as the only unknown: 

-v(S1, S2) 

P(S2)-P(S1)=X 2 (33) 



This new equation can be converted into a quadratic in X 2 and solved. The erroneous value of X 2 is 
then discarded and the correct value is substituted into Equation (28) to determine Z 2 . The value of Z 2 is 
then substituted into Equation (27) to determine Y 2 The sign of the dot product then reduces the quadrant 
ambiguity from 8 to 4. The 4 quadrant ambiguity is eliminated by specifying particular parameters as 
45 previously discussed. Figure 11 illustrates a flow diagram for the computations involved in the 
implementation of a two-state power and dot product solution for position. Figure 11 also goes on to 
illustrate the flow diagram for the computations involved in calculating orientation from the two-state 
power and dot product solution for position. The mathematical operations for calculating orientation are 
presented later. 



3. Three-state cross-product solution: 

The amplitude of the cross-product of two vectors is invariant under the orthogonal sensor orientation 

rotations and can therefore be used to determine position independent of sensor orientation. The 
65 orientation of the cross product is, in an absolute sense, also invariant under sensor orientation. However, 

the cross-product is referenced to the same coordinate frame as are the vectors used to generate it. 

Therefore, the cross-product of two sensor responses is referenced to the sensor coordinate frame. Since 

the orientation of the sensor is not known at this point in the signal processing, little use can be made of the 

orientations of the cross products. 
60 Inspection of the sensor-position-frame fields for the three axis excitation states equations (10), (1 1 ), 

and (12) shows that all have zero X components, which implies that both the Y and the Z components of the 

cross products of these vectors are zero. The resultant sensor-position-frame cross-product vectors are 

then 
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h (SI, S2)=& (S1)x& <S2) 
'(-sin a) (sin a sin (3)-(cos a) (cos a sin P) 

0 
0 

-sin p 
0 

. 0 



(34) 



(35) 



75 



j§3 (S2, S3)=i 3 (S2)xi3 (S3)=- 



C 2 



20 



and 



& (S3, SU-fe (S3)xj3 (S1)=- 



25 



cos a cos p" 
0 
0 

sin a cos 3 
0 
0 



(36) 



(37) 



30 



These three cross products will be rotated by yet unknown sensor orientation angles. However, the 
magnitudes of the cross-products (or the squared magnitudes) are unchanged, hence we can obtain, 
independent of sensor orientation, 



E(S1, S2)= |£(S1, S2) | 2 = 



C 4 

2 - sin 2 p 



35 



i(S2, S3)= |£(S2, S3) | 2 = cos 2 a cos 2 p 



40 



and 



P 4 



C 4 



(38) 



(39) 



E(S3, S1)= |£(S3, Si) I 2 = sin 2 a cos 2 p 



(40) 



45 



The solution for p, a 1f and p is in general similar to the solutions by other methods. First: 

p-C[=(Sl, S2)+E(S2, S3)+E(S3, SI)]" 1 ' 4 (41) 
so Substitution of p into (38) then produces: 

f> 2 



I P |=arcsin[ VE(S1, S2)] 



(42) 



55 



60 



and the ratio of (40) to (39) gives ambiguous azimuth: 

12 



a A =arctan[ 



E(S3, S1) 
E{S2, S3) 



(43) 



This position solution contains an eight-quadrant ambiguity than can be reduced to a two quadrant 
55 ambiguity through the signs of the dot products as shown in Table 3. 
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4. Two-state power and cross-product solution 

It is evident that P(S1) # P(S2), and =(S1, S2) produce three equations in the three unknown position 
parameters. To find position, first rearrange Equation (38) to produce 



P 4 

sin 2 P= E(S1, S2) 



(44) 



Substitution of this into the sum of Equations (14) and (15) leaves only p as an unknown: 



C 2 



P(S1)+P(S2)= 



={S1, S2) 



(45) 
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After p has been determined, | P | and a A can be determined in turn by substitutions into Equations (44) 
and (19), respectively. Quadrant ambiguity can be reduced from 8 to 8 by use of the sign of v(S1, S2). 

Orientation determination: 

Sensor orientation can be determined in a non-iterative manner from any two sensor output vectors 
corresponding to fields aligned with the source axes. These sensor output vectors are synthesized from the 
true sensor output vectors. One advantage of non-iterative orientation determination over iterative 
orientation determination is an increase in processing speed. Also, non-iterative orientation determination 
techniques are free from 'latch up' and allow a reduction in software complexity. 

The orientation rotations that convert the output of an equivalent zero-orientation sensor into the 
output of the true sensor 100 can be combined into a single matrix A, which can be expanded by using 
Table 1 to produce: 

A=I* Is I* 

cos 8 cos 4j cos 0 sin 41 -sin 6 

-cos <p sin i|i cos <p cos iJj 

+sin <|> sin 9 cos iJj +sin <p sin 6 sin \\t sin $ cos 9 (46) 

sin (p sin *(i —sin ip cos i|j cos <p cos 6 

+cos <p sin 9 cos i|j +cos <p sin 0 sin tjj 
Suppose that a source excitation products a response 

L(X)=[ 1, 0, 0] T 

in a zero-orientation sensor (i.e. the field at the sensor location has an X t -axis orientation. The output from 
the real sensor 100 is then 

i(X)=Ai 4 (X)=[a 11 , a 21 , a 3 il T , 

which is the first column of A. Similarly, the second and third columns of A represent the sensor responses 
to fields of YT-Z^axis orientations, respectively. 

If the normalized sensor output vectors corresponding to fields of X, -, Y r , and 2, -axis orientations can 
be synthesized, the elements of A will be known, and the angles 0, and <p can be determined. For 
example, using the sensor X-axis response to a field of Z r axis orientation, 



so 



0=-arcsin fi Jt (Z)=-arcsin a 13 



(47) 



The angles and <J> can now be determined by using the just-determihed value of 0 to cancel the sin 0 
and cos 6 factors in the responses corresponding to a t1 , a 12 , a 2 3i and a 33 . 

Errors in the range estimate and variation in the source power produce a multiplicative error common 
to all sensor output vectors. The effects of those errors can be avoided by determining orientation from 
ratios of sensor responses. Thus: 



i|j=arctan 



r f 6 , <Y) 1 

I =arctan 

L f s , <X) J 




(48) 
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to 



15 



20 



25 



30 



<(>=arctan 



f 5 y (Z) 



fsz (Z) 



=arctan 



a 2 3 



a 33 



(49) 



(Note that a four-quadrant inverse tangent will place ip and $ in the proper quadrant). Elevation 9 can 
be determined from 



6=arctan 



-f Sx (Z) 1 r -a 13 -| 

— =arctan 

_ f 5x (X)/cos 4» J L an/cos «jj J 



(50) 



or three similar ratios using a 13 , and a 12 , a 23 , or a 33 . A linear combination of all four ratios can also be used 
to minimize the effects of noise. 

While orientation is most simply determined using elements from all three synthetic sensor output 
vectors, inspection of the matrix A in Equation (46) shows that the information contained in any two 
columns is sufficient to determine all three orientation angles. Some flexibility is therefore possible in a 
noisy environment; e.g., orientation can be estimated from the two output vectors with minimum 
estimated noise. Alternatively, information from all three output vectors might be combined by linearizing 
the elements of A about the initial orientation estimates from Equations (48) r (49), and (50). Minimum- 
variance linear combinations would then be formed to improve the initial estimates. 

Near-field conditions, three-state excitation: 

Fields of X r , Y r , and Z r axis orientations are produced at the sensor location only when the sensor 100 
is located on the X v Y„ or Z, axis. The source excitation pattern is fixed to allow multiple sensors to derive 
position and orientation information from the same signals. In a near-field system with a three-state 
source-excitation pattern (United States Patent No. 4,054,881), the responses of the sensor to fields of X 1# 
Y„ and Z t -axis orientations can be synthesized from the real sensor responses, which span three- 
dimensional vector space. 

The three true sensor output vectors can be assembled into a 3x3 matrix Fg, which can then be written 

as: 



35 



fc(si) |i(S2) i.yss)]^ a Ql 

p 3 



A Q Ifa(Sl) ii(S2) : i,(S3)] 



(51) 



40 From the above equation, it is evident that the desired matrix A of synthetic responses can be obtained 
as: 



45 



A=- 



(52) 



The coupling matrix 6" 1 is computed using estimated values a and (3, Note that actual matrix inversion 
is unnecessary, since 



so 
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where 



Q-^d-a T_ p § Tp JJ-^I-o I- P S" 1 la la 
1 0 0 1 



(53) 



0-2 0 
L 0 0-2 



(54) 



Far-field conditions, three- or two-state excitation: 

While all orientation information is contained in any two sensor output vectors, synthesis of the 
desired sensor responses {i.e., the matrix A) requires a three-dimensional set of basis vectors. However, 
under fa r-fiel decoupling conditions, the coupling matrix S (Equation (3)) is degenerate (rank 2). Therefore, 
the inverses S \ hence Q" 1 , (Equations (53) and (54)) do not exist and Equation (52) cannot be used directly 
to synthesize orientation matrix 

The orientation matrix A can, however, be synthesized by using the cross product of two non-colinear 
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sensor output vectors to provide the necessary third linearly independent vector. Suppose that the 
responses to state S1 and state S2 excitation are available and since orthogonal rotations preserve the 
angles between vectors, 

5 i 5 (CP)=i 5 (S1)xi s (S2) 

=[Ai,(S1)]x[Ai(S2)]=A[t(S1)xi(S2)] (55) 

The cross products can be used in place of the vectors produced by the third excitation state in forming 
;o the matrices fs and JF*, thus 

Fs-lfglSD i(S2) i5(CP)]=A&<S1) £<S2> £<CP)]=A £, (56) 

The vectors f* (S1) and j4(S2) and hence f^CP) can be calculated from the estimated position. The 
is vector fg(CP) can be calculated from the two sensor output vectors. The estimated orientation matrix A can 
then be determined from: 

A=£s JET 1 (57) 

20 Orientation angles are then determined as discussed previously. For computational convenience, it 

may be desirable to multiply ifSD, fs(S2), j^SU, and.£,(S2) by p/C before computing the cross products; 

this results in cross-product vectors and sensor output vectors having roughly the same magnitude. 
Matrix inversion can be avoided by an alternative method of determining orientation. Linear 

combinations of the two field vectors at the sensor location form two orthogonal synthetic field vectors; 
25 application of the same coefficients to the two sensor output vectors produces the analogous synthetic 

response vectors. 

The two synthetic sensor response vectors are placed in the second and third columns of a matrix A', 
and orientation angles t|i', 9', and <p' are determined from those two columns. Those orientation angles are 
defined in the same manner as are qj, 8, and but are referenced to the Y'-Z' coordinate frame formed by 
30 the two synthetic field vectors. The orientation (a', 0\ 0') of the Y'-Z' axes with respect to the X^Y^Z, 
axes is then determined by multiplying direction cosine vectors. The matrix A whose orientation angles are 
referenced to the X, -Yt -Zt coordinate frame is then: 

A=I/ Te' X/ 1/ TV L' (58) 

35 

This method may allow a reduction of computation time in some applications. 

If the sensor lies in the X, -Y t plane, two sensor outputs will be colinear, and orientation cannot be 
determined from fs<S1) and ig(S2) alone. This suggests that for the most general allowable orientations, 
three-state operation should be used so that there will always be some basis against which orientation can 
40 be determined. 

In far-field operation, the cross-product of any two field vectors at the sensor location is oriented 
radially away from the source, see Equations (35), (36), and (37). Transmission of such a vector is physically 
impossible, hence the cross-product is not a direct substitute for the third-axis (S3) excitation. 

The cross-product is nonetheless useful in determining large-angle orientation for far-field operation. 

45 Linear combinations of the two real sensor responses and their cross product can produce synthetic 
responses to non-physically realizable source-frame Y n - and Z, -oriented fields. The synthetic responses 
25(S2) andas(S3) thus determined then yield source-frame orientation angles when used in the large-angle 
orientation algorithms. The coefficients required to implement the two linear transformations are elements 
of the inverse of a matrix composed of the two computed field vectors [£j(S1) and.&(S2)] and their cross 

so product f^CP). 

It should be noted that the equations which have been derived in the present disclosure are based 
upon electromagnetic fields generated by electrically short (0.1 A or shorter) dipoles. Longer dipoles and 
arrays have different field patterns that will make these equations invalid. If such antenna are to be used, 
appropriate equations should be derived. 
55 it will be evident to an expert in the art how to program the analysing means 50 to apply the methods 
and the equations listed above in any particular case. 

Signal formats: 

The signal format used by the transmitter must be designed to allow the user to determine his position 
so and orientation. The geometric computations discussed in the previous section establish some 
informational parameters which are advantageously met by the format. First, it must allow airborne 
components 20 to determine the amplitude of the signal induced in a given receiving antenna (22, 23 or 24) 
by each transmitting antenna (41, 42 or 43). Secondly, it advantageously provides one-way data 
transmission capability so airborne components 20 will know the power, (i.e. strengths), of the transmitted 
65 signal. Thirdly, to facilitate communication of both the previous information parameters, the signal may 
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include a timing reference, and all signal components should advantageously be coherent with this 
reference. The timing signal is used to characterize the polarity of the transmitted signal. If this timing 
signal is omitted, there is an increase in the ambiguity as to position and orientation of the remote object. 
Of course, independent information sources may be used to remove this ambiguity. For example, 

5 navigation aids may be used to determine the quadrant (i.e. northeast, southeast, northwest or southwest) 
of the remote object with respect to the radiating means; altimeters may be used to determine the relative 
height of the remote object with respect to the radiating means, which may be located at the top of the hill. 

There are endless varieties of formats which can meet the above requirements. However, it is 
additionally desirable that the signal format allow easy acquisition by the user as he approaches the 

to landing point. Simplicity in the receiver is also quite desirable. Four possibilities for transmitted signal 
format are: 

1. Frequency division multiplexing (FDM) 

In this format, each transmitting antenna in the array is assigned a particular, different frequency. 
rs Measurements of the information parameters can be simply the outputs of integrators corresponding to 
the three frequencies. The carriers are of constant phase and therefore easily acquired by a phase-locked 
loop with an appropriate time constant. 



20 2. Time division multiplexing (TDM) 

In TDM, only one dipole of the ground antenna array is excited at a time. Simplification of transmitter 
and receiver is possible because circuits can be time-shared. However, data transmission is more complex, 
and moving aircraft must interpolate between measurements to provide the equivalent of simultaneous 
measurements. This is an advantageous mode of multiplexing for the location of a relatively slow moving 

25 vehicle. 



3. Phase division multiplexing (nutation) 

Proper excitation of the three dipoles will produce the equivalent of physical nutation of a single dipole 
30 antenna. This may be accomplished by excitation of the 2 dipole with an unmodulated carrier and the X 
and Y dipoles by carriers with high frequency amplitude modulation by sine and cosine waves, 
respectively. While this results in a beacon-like signal, that property is not actually utilized in the position 
and orientation calculations. What is utilized is that the radiated electromagnetic field from each of the 
ground dipole antennas can be identified by either an unmodulated carrier or a phase difference between 
35 the modulation envelope of two modulated carriers. If a nutating field is used there is no requirement that 
the axis of nutation be positioned along a line between the radiating and the receiving antennas. The 
position and orientation of the remote object can be determined regardless of the orientation of the axis of 
nutation. 

40 

4. Spread spectrum multiplexing 

To accomplish spread spectrum multiplexing, each transmitted signal is assigned a unique code 
sequence which shifts carrier frequency, carrier phase (or both). Reception is accomplished by using 
identical code sequences to remove the modulation. The codes assigned to the three antennas are 

45 designed not to cross correlate and thus make possible measurement of individual signals. However, 
acquisition is typically more difficult, both because of the absence of a carrier component and because the 
code timing must be acquired, in addition to the carrier. 

The selection of carrier frequency, spread-spectrum chipping rate, data frequency, transmitter power, 
antenna size, and other parameters is naturally dependent upon this application. Carrier frequencies in the 

so 200 to 3000 MHz range will generally be found suitable. For these carrier frequencies, chipping rates of 100 
kHz to 10 MHz are practical. Consequently, data frequencies of 10 kb/s to 1 Mb/s are possible. If 
frequency-division multiplex is used, carrier frequencies should be separated by 10 kHz to 100 kHz to allow 
for Doppier shifts. 

An analogous technique can be used to determine the position and orientation of a two-axis sensor 
55 with respect to a three-axis source. Use of the equations given here produces source position and 
orientation relative to the sensor. A simple series of transformations converts this information to the 
desired coordinate frame. 

While the present invention has been described in terms of a long distance landing aid system, it can 
also be useful in applications such as airdrop guidance and control, collision avoidance, target handoff, and 
60 refueling and station keeping. 

Various modifications and variations will no doubt occur to those skilled in the various arts to which 
this invention pertains. As discussed above, the signal format may be chosen from any of numerous 
alternatives. Additionally, the particular parameters of the transmitting and receiving apparatus will 
depend upon the particular application. Systems for either longer or shorter ranges can be designed by 
55 appropriate choices of parameters. 
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Claims 

1. Apparatus for determining the position and orientation of a remote object relative to a reference 
coordinate frame comprising: a plurality of radiating means having orthogonal components (41, 42, 43) 

5 centered about the origin of the reference coordinate frame; means (30) for applying to the plurality of 
radiating means electrical signals which generate a plurality of electromagnetic fields, the electromagnetic 
fields being distinguishable from one another; a plurality of receiving means disposed on said remote 
object, said receiving means having orthogonal components (22, 23, 24) for detecting and measuring 
components of the electromagnetic fields; and analyzing means (50) associated with the receiving means 

70 for converting the received components of the electromagnetic fields into remote object position and 
orientation relative to the reference coordinate frame with at most two ambiguous combinations of 
orientation and position, the analyzing means operating open loop with respect to the radiating means; 
characterized in that one of the said pluralities consists of only two orthogonal components and the other 
consists of three components; that the apparatus is adapted for operation in a region in which the far-field 

75 components of the electromagnetic fields are substantially greater in magnitude than the near-field 
components; and that the orientation is determined by a method which includes obtaining from the 
received components two of three vectors representative of the orientation of the remote object and 
deriving the third vector by forming the cross product of the two known vectors. 

2. Apparatus as claimed in Claim 1 in which each of the said radiating means and the said receiving 
20 means is a mutually orthogonal, electric-dipole antenna. 

3. Apparatus as claimed in Claim 1 or Claim 2 in which the electromagnetic fields are distinguishable by 
means of time, frequency, phase or spread spectrum multiplexing of the said electrical signals. 

4. Apparatus as claimed in any of Claims 1 to 3 in which the electrical signals further include 
information relating to at least one of: transmitter identification, electromagnetic field distortion 

25 corrections, locations of obstacles, location of a landing site relative to the radiating means and wind 
direction. 

5. Apparatus as claimed in any of Claims 1 to 4, which includes a monitor station at a fixed location 
spaced from the radiating means and having receiving means for detecting the radiated electromagnetic 
fields and being coupled to the transmitting means for providing feedback to the transmitting means 

30 characterizing the transmitted fields. 

Patentanspruche 

1. Einrichtung zum Bestimmen der Position und der Richtung eines fernen Objektes relativ zu einem 
35 Bezugskoordinatensystem mit: Emer Mehrzahl von Abstrahlungsanordnungen mit orthogonalen 

Komponenten (41, 42, 43), die etwa in dem Ursprung des Bezugskoordinatenssystems angeordnet sind; 
einer Anordnung (30) zum Zufuhren von elektrischen Signalen, welche eine Mehrzahl von elektrischen 
Feldern erzeugen an die Mehrzahl der Abstrahlungsanordnungen, wobei die elektrischen Felder 
voneinander unterscheidbar sind; einer Mehrzahl von an dem fernen Objekt angeordneten Empfangs- 

40 anordnungen, die orthogonale Komponenten (22, 23, 24) zum Erkennen und Messen von Komponenten 
der elektrischen Felder aufweisen; und einer der Empfangsanordnung zugeordneten Ana lysierein richtung 
(50) zum Umsetzen der empfangenden Komponenten der elektrischen Felder in die Position und die 
Richtung des fernen Objektes relativ zum Bezugskoordinatensystem mit hochstens zwei zweifelhaften 
Kombinationen der Richtung und der Position, wobei die Anaiysieranordnung bezuglich der Abstrahlungs- 

45 anordnung entsprechend einer offenen Schleife arbeitet, dadurch gekennzeichnet, daft eine der Mehr- 
zahlen nur aus zwei orthogonalen Komponenten besteht und die andere aus drei Komponenten besteht, 
daft die Einrichtung fur ein Arbeiten in einem Gebiet angepaftt ist, in dem die Fernbereich-Komponenten 
der elektromagnetischen Felder hinsichtlich ihrer Amplituden wesentlich grofter sind als die Nah- 
bereich-Komponenten und daft die Richtung durch ein Verfahren bestimmt wird, welche das Erhalten von 

so zwei oder drei Vektoren, die die Richtung des fernen Objektes darstellen aus den empfangenen 
Komponenten und die Ableitung des dritten Vektors durch die Bildung der Kreuzmultiplikation der beiden 
bekannten Vektoren umfaftt. 

2. Einrichtung nach Anspruch 1, dadurch gekennzeichnet, daft jede der Abstrahlungsanordnungen und 
jede der Empfangsanordnungen als eine wechselseitig orthogonale elektrische Dipolantenne ausgebildet 

55 ist. 

3. Einrichtung nach Anspruch 1 oder 2, dadurch gekennzeichnet, daft die elektrischen Felder 
voneinander durch die Zeit, Frequenz, Phase oder das Multiplexen eines gespreizten Spektrums der 
elektrischen Signale unterschiedbar sind. 

4. Einrichtung nach einem der Anspruche 1 bis 3, dadurch gekennzeichnet daft die elektrischen 
60 Signale weiterhin eine Information enthalten, die sich auf wenigstens eines der folgenden bezieht: 

Senderidentifikation, Korrekturen oder Verzerrung des elektrischen Feldes, Lagen von Hindernissen, Lage 
eines Landeplatzes relativ zur Abstrahlungsanordnung und die Windrichtung. 

5. Einrichtung nach einem der Anspruche 1 bis 4, dadurch gekennzeichnet, daft sie eine Anzeigestation 
enthalt, die an einem festen, von der Abstrahlungsanordnung entfernten Ort angeordnet ist und die eine 

65 Empfangsanordnung zum Erkennen der abgestrahlten elektromagnetischen Felder aufweist und die mit 
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der Sendeanordnung verbunden ist, urn eine die ausgesendeten Felder kennzeichnende Ruckfuhrung an 
die Sendeanordnung durchzufuhren. 

Revendications 

1. Appareillage pour determiner la position et I'orientation d'un objet eloigne par rapport a une 
systdme de coordonnees de reference comprenant: une pluralite de moyens de rayon nement comportant 
des composants orthogonaux (41, 42, 43) centres sur I'origine du systeme de coordonnees de reference; 
des moyens (30) pour appliquer a la pluralite de moyens de rayonnement des signaux electriques qui 
engendrent une pluralite de champs electromagnetiques, les champs 6lectromagnetiques pouvant etre 
distingues I'un de I 'autre; une pluralite de moyens recepteurs disposes sur ledit objet eloigne, lesdits 
moyens recepteurs comportant des composants orthogonaux (22, 23, 24) pour detecter et mesurer des 
composantes des champs electromagnetiques; et des moyens analyseurs (50) associes aux moyens 
recepteurs pour convertir les composantes de champs electromagnetiques regues en une position et une 
orientation de I'objet eloigne par rapport au systeme de coordonnees de reference, avec au plus deux 
combinaisons ambigues d'orientation et de position, les moyens analyseurs operant en boucle ouverte par 
rapport aux moyens de rayonnement; caracterise en ce qu'une desdites piuralites comprend seulement 
deux composants orthogonaux et I'autre comprend trois composants; en ce que I'appareillage est adapts 
pour operer dans une region ou les composantes de champ eloigne des champs electromagnetiques ont 
une grandeur sensiblement plus grande que les composantes de champ proche; et en ce que I'orientation 
est determinee par un procede qui consiste a obtenir, a partir des composantes regues, deux des trois 
vecteurs representant I'orientation de I'objet 6loign6 et a d6duire le troisieme vecteur par formation du 
produit vectoriel des deux vecteurs connus. 

2. Appareillage tel que revendiquS dans la revendication 1, dans lequei chacun desdits moyens de 
rayonnement et desdits moyens recepteurs est une antenne a dipoles electriques mutuellement 
orthogonaux. 

3. Appareillage tel que revendique dans la revendication 1 ou la revendication 2, dans lequei les 
champs electromagnetiques sont distinguables au moyen d'un multiplexage en temps, en frequence, en 
phase ou en spectre d'etalement desdits signaux electriques. 

4. Appareillage tel que revendique dans Tune quelconque des revendications 1 a 3, dans lequei les 
signaux electriques contiennent en outre une information concernant au moins un des paramdtres 
suivants: identification d'emetteur, corrections de distorsion de champ electromagnetique, emplacements 
d'obstacles, emplacement d'un site d'atterrissage par rapport aux moyens de rayonnement et a la direction 
du vent 

5. Appareillage tel que revendique dans Tune quelconque des revendications 1 a 4, qui comprend un 
poste de contrdle en un endroit fixi et espac§ des moyens de rayonnement, et comportant des moyens 
recepteurs pour detecter les champs electromagnetiques rayonnes, et qui est relie aux moyens emetteurs 
pour fournir a ces moyens emetteurs une retroaction caracterisant les champs emis. 
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